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Coupled map lattices of non-hyperbolic local maps arise naturally in many physical situations 
described by discretised reaction diffusion equations or discretised scalar field theories. As a proto- 
type for these types of lattice dynamical systems we study diffusively coupled Tchebyscheff maps 
of 7V-th order which exhibit strongest possible chaotic behaviour for small coupling constants a. 
We prove that the expectations of arbitrary observables scale with in the low-coupling limit, 
contrasting the hyperbolic case which is known to scale with a. Moreover we prove that there 
are log-periodic oscillations of period log A''^ modulating the y^-dependence of a given expectation 
value. We develop a general fst order perturbation theory to analytically calculate the invariant 
1-point density, show that the density exhibits log-periodic oscillations in phase space, and obtain 
excellent agreement with numerical results. 
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I. INTRODUCTION 



Coupled map lattices (CMLs) as introduced by Kaneko 
and Kapral [l[ [2] are a paradigm of higher-dimensional 
dynamical systems exhibiting spatio-temporal chaotic 
behaviour. There is a variety of applications for CMLs 
to model hydrodynamical flows, turbulence, chemical re- 
actions, biological systems, and quantum field theories 
(see e.g. reviews in [3, The analysis of chaotic CMLs 
is often restricted to numerical investigations and only a 
few analytical results are known. A notable exception is 
the case of hyperbolic maps (maps for which the absolute 
value of the slope of the local maps is always larger than 
1) for small couplingo. In this case a variety of analytical 
results exists 0, H, U 111 that guarantee the existence of 
a smooth invariant density and ergodic behaviour. The 
situation is much more complicated for non-hyperbolic 
maps which correspond to the generic case of physical in- 
terest. Here much less is known analytically, though some 
promising steps have been made 0, HO, El, [H, [13, ■ 

In this paper, we present for the first time analytical 
results corresponding to a non-hyperbolic situation and 
calculate invariant densities explicitely for the case of lo- 
cally fully developed chaos that is diffusively coupled. 
We study local maps given by iV-th order Tchebyscheff 
polynomials. In the uncoupled case these maps are con- 
jugated to a Bernoulli shift of N symbols. In the cou- 
pled case, this conjugacy is destroyed and the convential 
treatments for hyperbolic maps does not apply, since the 
Tchebyscheff maps have — 1 critical points where the 
slope vanishes, thus corresponding to a non-hyperbolic 
situation. From the physical point of view, the non- 
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hyperbolic case is the most interesting one. For exam- 
ple, it has been shown that these types of non-hyperbolic 
CMLs naturally arise from stochastically quantised scalar 
field theories in the anti-integrable limit They can 
serve as useful models for vacuum fluctuations and dark 
energy in the universe [isj . Other applications include 
chemical kinetics as described by discretised reaction dif- 
fusion dynamics 

The case of two coupled Tchebyscheff maps was previ- 
ously studied in (l6| using periodic orbit theory. Here we 
consider infinitely many diffusively coupled Tchebyscheff 
maps on one-dimensional lattices with periodic boundary 
conditions, and apply Perron- Frobenius and convolution 
operator techniques. Our analytical techniques will yield 
explicit perturbative expressions for the invariant 1-point 
density for small couplings a. We will prove that the 
density exhibits log-periodic oscillations of period log A^^ 
near the edges of the interval. Our explicit result for the 
invariant density will allow us to calculate expectations 
of arbitrary local observables. We will prove that expec- 
tations of typical observables scale with (rather than 
with a as for hyperbolic coupled maps). We also show 
that there are log-periodic modulations of expectation 
values when the parameter a is changed. Our results 
seem to be typical for local maps with one or several 
quadratic maxima which are in a fully developed chaotic 
state. Other types of maps may of course generate dif- 
ferent types of behavior for a ^ [ll[ . 

This paper is organised as follows: In Sec. II we in- 
troduce the relevant class of coupled map lattices and 
briefly explain their physical relevance. In Sec. HI we 
give some numerical results for the scaling behaviour of 
these non-hyperbolic systems. In Sec. IV we present our 
analytical results for the invariant density which we use 
to rigorously prove the scaling behaviour. 
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II. THE CLASS OF SYSTEMS 

Consider a scalar field (p(af, t) described by an equation 
of the form 



d_ 

dt 



if = D ■ Aip + V'{ip). 



(1) 



Here D is a diffusion constant and V is some suitable 
potential. These types of equations occur in many dif- 
ferent areas of physics. They describe reaction diffu- 
sion systems, Ginzburg-Landau type of models, nonlinear 
Schrodinger equations, stochastically quantised field the- 
ories, etc. In the one-dimensional case the Laplacian A 
is just given by d'^/dx^. In many cases there is a fun- 
damental length and time scale below which the above 
continuum theory is not valid anymore. For example, for 
stochastically quantised field theories this is the Planck 
length Writing t = nr, x = i6 where n and i are inte- 
gers, and fixjt) = Pmax^li where is a dimensionless 
field and Pmax is some constant with the same dimension 
as the field, the discretised Eq. ([T]) can be written in the 
form 



where the local map T is given by 
T($) = $ 



(1 - a)pn 



-V^'bmax*), 



(2) 



(3) 



and the dimensionless coupling a is given by a = 2Dt/S^. 
For various applications of these types of coupled map 
lattices, see The important point is that generically 

these types of models can lead to non-hyperbolic local 
maps T. For example, a i^s'^-theory described by a dou- 
ble well potential V{ip) with a sufRciently strong quartic 
term leads to cubic maps with two inflection points where 
the slope T'(<I>) vanishes. It is thus important to under- 
stand the generic behaviour of non-hyperbolic coupled 
map lattices. 

We are particularly interested in cases where the local 
map exhibits strongest possible chaotic behaviour. The 
negative 3rd order Tchebyscheff map 



n+l 



-4$3 + 3$„ 



(4) 



on the interval $ G [—1,1] is such an example. It is 
conjugated to a Bernoulli shift of 3 symbols and can be 
obtained in our context from the potential 



V{ip) 



1 



2 



const. 



(5) 



In a similar way, one can construct potentials that lead 
to positive and negative Tchebyscheff maps of arbitrary 
order N H. 



III. OBSERVED SCALING BEHAVIOUR 

Let us now study CMLs of type ^ where T = Tjy is 
a Tchebyscheff map of order N . One observes nontrivial 
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FIG. 1: The function (^'(loga) of Eq. (0 witli = 
Viv($) as function of a for A'^ = 2, 3 in a double-logarithmic 
plot with basis 7V^. 



scaling behaviour for small values of the coupling a that 
is significantly different from that of hyperbolic systems. 
Consider an arbitrary test function of the local it- 
erates. Assuming ergodicity, expectations (/i($))a for a 
given parameter a are numerically calculable as time av- 
erages 



M J 



(6) 



For a 



one numerically observes the scaling behaviour 
{h{^))a~{hm^ = V^-F^''\\oga) (7) 



where is a periodic function of log a with period 

logA^^. Examples are shown in FIG. [TJ The choice 
h{^) = V±jv(<i)) := t/ r7v(<i>)d$ is important in the 
physical applications to estimate the vacuum energy gen- 
erated by the chaotic field theory under consideration [41 . 
Generally the function F — F\h] \s a. functional of the 
chosen test function h. 

The above log-periodic scaling is observed for arbitrary 
test functions h and hence is a general property of the in- 
variant 1-point density Pa{^) of the CML for given small 
couplings a. In fact one observes that there is not only 
scaling in the parameter space a but also in the phase 
space $. This is shown in FIG.H] 

Near the left edge of the interval [—1,1] one may write 
$ = ay — 1 and observe the scaling behaviour 



Pa{ay - 1) = a ^/^.g(y) 



(8) 



where the function g is independent of a for small a. At 
the right edge, writing <I> = 1 — ax, one observes 

Pa{l - ax) = po(l - ax) -f ]^a~^'^x~^f{x) (9) 

where / is independent of a for small a. Moreover, / 
exhibits log-periodic oscillations 



J{N^x) = fix) 



(10) 
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FIG. 2: Log-periodic oscillations of the rescaled invariant den- 
sity f{x) as observed for N — 2 (top) and TV = 3 (bottom) 
and different values for a. 



over a large region of the phase space (see FIG. [2]). The 
number of oscillations is approximately — log^2 (a). 

For N odd the density is symmetric and the behaviour 
at the left and right edge is the same, whereas for TV even 
there is no left-right symmetry. 

In the following section we will analytically prove the 
above numerical observations and provide explicit formu- 
las for the functions / and g for a given local Tchebyscheff 
map T/v. 



IV. PERTURBATIVE CALCULATION OF THE 
INVARIANT DENSITY 

Our method is based on a perturbative treatment of 
the Perron-Frobenius operator of a perturbed local map 
and on convolution techniques. In a first approximation, 
the neighbouring lattice sites can be regarded as pro- 
ducing independent random noise with density poi^) — 
(1 — $^)^^/^/7r. While this picture works well in the mid- 
dle of the interval [—1,1], in the vicinity of the edges ±1 
one has to take into account nontrivial nearest neighbour 
correlations, due to the non-hyperbolicity of the map. A 
first-order approximation of the density can then be fur- 
ther iterated to yield results of better precision. In each 
step we integrate over 2-point functions describing the 
joint probability of neighboured lattice sites to obtain 



the marginal distribution at a single lattice site. A de- 
tailed description of our calculations is out of the scope 
of a short letter, they will be pubHshed elsewhere [rn |. 
Our final result is that for N = 2 one obtains in leading 
order of ^/a at the left edge Paiay '''^ ~ "'"^^ 
where 



1) 



1 



1) « p^iay 
Pa{4>+)d(f>+Po{(f)-)d(f)^ 



1) 



(11) 



/2a 7 y/y-l + {<l>++cb-)/2 
At the right edge one obtains by iterating our scheme q 



times pa(l 



ax) 



Pa 



1 



-.1 Pa 



4P7rV2a 



(1 — ax) where 

Pa{<P+)d(f>+PQ{(l)-)d(j)- 

Here the function r2((/)) is defined as follows: 

^2,(0) - 1 



1 P 

2 ^ 

q=0 



22q 



■-) 

(12) 



(13) 



The limits of the two integrations in Eqs. pT|) and (fT^ 
are given by the condition that \4>±\ < 1 and that the 
argument of the square root should always be positive. 

A simple way to numerically evaluate our formulas is 
to replace the double integrals by ergodic averages of iter- 
ates of the uncoupled Tchebyscheff map. So for example 
we may evaluate the density at the left edge as 



P^°\ay- 



M 



0(2/ - 1 -K$„^ + J/2) 



(14) 



and similar formulas apply to the right edge. 
Eq. pT|) can also be written as 



p^^Hay-i) = 



Poo{z) 



1 



2 



/2a 



PQo{z)dz 

- 1 + 2 



(15) 



KiVl~z^)d{l~-z^) (16) 



where K{x) is the complete elliptic integral of the first 
kind. From the above formula it is obvious that the den- 
sity is not differentiable at y = 1 and y = 2 for arbi- 
trarily small a (compare FIG. [3]). By iteration of the 
Perron-Frobenius operator one can show that these two 
non-analytic points generate an entire cascade of such 
points at the right edge of the interval. 

Note that the sum over the terms in Eq. (|12p con- 
verges rapidly with increasing q. In practice, a few terms 
are sufficient to obtain perfect agreement with the nu- 
merical histograms. FIG. [3] shows how well our analyt- 
ical results (|lip and agree with the numerics. The 
agreement is so good that the analytic curves (given by 
dashed-dotted lines in FIG. [3|) are not visible behind the 
data points. 

Along similar lines, we obtain for TV = 3 

(p)/. ^ 2 f po{(l)+)d(l)+po{(f>^)d(f>^ 

Pa'{^~ax'- 



9PST:V2aJ ^/x/9p + {(/)+) + 



(17) 
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FIG. 3: Scaling behaviour of the invariant 1-point density 
near the left and right edge of the interval [-1,1] for N = 2 
and a — 0.00007625. Shown is the density at the left (upper 
diagram) and right edge (lower diagram) in comparison with 
the exact results (dashed-dotted curve, covered by the nu- 
merical result) and the density po{4') (dotted curve) together 
with several thresholds, i.e. points with nonanalytic behavior 
(dashed lines). 



with 



9=0 



91 



The density is symmetric, i.e. 

For general N it is natural to conjecture that 



(18) 



(19) 



Pai(t>+)d(t)+Poi(j)-)d(l). 



' N 



1 

-T 

2 ^ 

9=0 



(21) 



(20) 



With the above equations we have derived explicit for- 
mulas for the scaling functions / and g in Eqs. ^ and 
It is easy to check from this integral representation that 
/ indeed satisfies property (fTU)) . Using Eqs. ([5|), (O 
and (|10p it is also possible to prove the general scaling 
relation ([7]) for arbitrary test functions h. Due to space 
restrictions we do not describe the details here but refer 
to a longer version [l7|. The function F^^^ is essentially 
given by a suitable integral of the observable h folded 
with /. 



V. CONCLUSION AND OUTLOOK 

We have derived in a perturbative way the 1-point den- 
sity for coupled Tchebyscheff maps and rigorously proved 
the existence of several interesting scaling phenomena 
that are numerically confirmed. Although our formulas 
were worked out for the explicit example of Tchebyscheff 
maps, many of our techniques and results are expected 
to hold in a similar way for other non-hyperbolic systems 
exhibiting strongly chaotic behaviour. Scaling with y/a 
and log-periodic oscillations seem to be typical phenom- 
ena for nonhyperbolic dynamical systems for which the 
local map has one or several quadratic maxima and is in 
a fully developed chaotic state. 
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